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Abstract 

The thermal averaged real-time propagator of a Dirac fermion in a static 
uniform magnetic field B is derived. At non-zero chemical potential and tem- 
perature we find explicitly the effective action for the magnetic field, which is 
shown to be closely related to the Helmholz free energy of a relativistic fermion 
gas, and it exhibits the expected de Haas - van Alphen oscillations. An effec- 
tive QED coupling constant at finite temperature and density is derived, and 
compared with renormalization group results. We discuss some astrophysical 
implications of our results. 
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1 Introduction 



Large magnetic fields B can be associated with certain compact astrophysical objects like 
supernovae |M [| where B = (9(10 10 )T, neutron stars ^J] where B = (9(10 8 )T, or white 
magnetic dwarfs 0, [IJ in which case S = (9(10 4 )T. (As a reference the electron mass in 
units of Tesla is m 2 e /e = 4.414- 10 9 T.) It has recently been argued that a plasma at thermal 
equilibrium can sustain fluctuations of the electromagnetic fields. In particular, for the 
primordial Big-Bang plasma the amplitude of magnetic field fluctuations at the time of the 
primordial nucleosynthesis can be as large as B = 0(1O 1O )T ||. Furthermore, a model for 
extragalactic gamma bursts in terms of mergers of massive binary stars suggests magnetic 
fields up to the order B = (9(10 13 )T 0. A more speculative system where even larger 
macroscopic magnetic fields can be contemplated are superconducting strings ||. Here 
one may conceive fields as large as B > O(10 14 )T. It has also recently been suggested 
that due to gradients in the Higgs field during the electroweak phase transition in the 
early universe very large magnetic fields, B = O(10 19 )T, may be generated @. If one 
encounters magnetic fields of this order of magnitude the complete electroweak model 
has to be considered and the concept of electroweak magnetism becomes important (for 
a recent account see e.g. Ref.JlO]). In the present paper we consider, however, magnetic 
fields such that calculations within QED are sufficient. A shorter version of this report 
has been published elsewhere [ |TT|] . 

In many of these systems one has to consider the effects of a thermal environment 
and a finite chemical potential. In this paper we derive the appropriate effective fermion 
propagator and the effective action in QED for a thermal environment treated exactly in 
the external constant magnetic field but with no virtual photons present, i.e. we consider 
the weak coupling limit. 

Calculations of the QED effective Lagrangian density in an external field have been 
attempted before either at finite temperature [|T2| , |I~3|j or at finite chemical potential Jl4| . 
In the latter case [II] the effective action is not complete but the correct form is presented 
here. At finite chemical potential and for sufficiently small temperatures, the QED ef- 
fective action should exhibit a certain periodic dependence of the external field, i.e. the 
well-known de Haas - van Alphen oscillations in condensed matter physics. This was 



not obtained in Ref.|L4|. Elsewhere, the radiative corrections to the anomalous magnetic 
moment has been estimated in the presence of large magnetic fields and it was argued 



that they are extremely small |T5| [16 
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By making use of the effective action we derive the effective QED coupling as a function 
of the external field, the chemical potential and temperature. In a future publication we 
will discuss the fermion self-energy and radiative corrections to the electrons anomalous 
magnetic moment, in terms of the formalism derived here fI7 |. 



2 Thermal propagators in the Furry picture 



We consider Dirac fermions in the presence of an external static field as described by the 
vector potential A^. Using static energy solutions we may represent the second quantized 
fermion field in the Furry picture Hl8| . It is given by 

M/(x,t)=y:6 AK ^: ) (x,t)+4^t ) (x,t) , (2.i) 

X.K 

where A is a polarization index, k denotes the energy and momentum (or other) quantum 
numbers (discrete and/or continuous) needed in order to completely characterize the 
solutions, and (±) denotes positive and negative energy solutions of the corresponding 
Dirac equation, 

(tP-m)^\x,t)=0 , (2.2) 

where = d^ + ieA^ is the covariant derivative. The creation and annihilation operators 
satisfy the canonical anti-commutation relations 

{d\> K >, 4J = <5a'a<5 k ' k = {b x > K >, b[ K } , (2.3) 

while other anti-commutators are zero. The completeness relation 

£ <IV, *)V>S(*, t) + < } a V, t)< } fe (x, t) = 5 ab 5 3 (x' - x) , (2.4) 

where ^ita denotes the a-component of the Dirac spinor if)^ , leads to the canonical 
ant i- commutation relations for the fields 

{^ a (x',t),^(x,t)} = ^ 3 (x'-x) . (2.5) 

In vacuum, the fermion propagator iSf{x'\ x\m) is defined by 

iS F (x';x\m) = (0|T (*(x, t')V(x,t)) |0) = 

0(*'-*)£v£V^itW) -0(*-O£4;V,^iW) , (2.6) 

\k Aft 
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where the conjugated spinor is given by ijj^ = (ipi^Ylo- Since ip^\x,t) satisfies 
the Dirac equation, only the time derivative acting on the step functions gives a non-zero 
contribution, so one finds that 



{ilp — v(i)Sf{x'\ x\m) = 1-5 (x' — x) 



(2.7) 



The real-time propagator at finite temperature T and chemical potential /a, denoted 
by {iSf{x'; x\rri)) can now be obtained by the following reasoning. Let fp{oj) denote 
the Fermi-Dirac thermodynamical distribution function 

1 



(2.8) 



exp(/3(iu -//)) + 1 

where f3 is the inverse temperature, fi the chemical potential and u is the energy of 
the quantum state under consideration. A particle can propagate forward in time in a 
state which is unoccupied by thermal particles, whereas a hole in the occupied states can 
propagate backwards in time . We can therefore write 



{iSf{x'; x\m))p ii 



0(t' -*)([!- /^)]V£V/)V>itW) + 



E 

X.K 



(+)/ 



We can now extract the vacuum part of the propagator Eq^^) and write 

(iSF(x'] x\m))fan = iSF^x'; x\m) + iSp{x'\ x\m) , 
where the thermal part iSp(x';x\m) is defined by 

Sp(x'',x\m) = 



(2.9) 



(2.10) 



and where we have defined the distribution 



X.i 



fr{E K ) = l-f+(-E K ) 



(2.11) 



(2.12) 



Notice that there is no time-ordering in Sp(x'; x\m) despite the fact that the time- 
ordering in Eq. ( |2.9|) is non-trivial. The thermal propagator Eq. fl2.10D therefore also triv- 
ially satisfies Eq. (|2.7|) . These considerations can, of course, easily be extended to treat 
particles with Bose-Einstein statistics as well. 
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The result Eq.( |2.11| ) can also be derived from an explicit calculation using the second- 
quantized field operators and appropriate thermal averages, i.e. we use Wicks theorem 



T(^(x',t')*( x ^)) =iS F (x';x\m) + : ^(x', t')^(x,t) : 
where the last term corresponds to a normal ordering. We then obtain 
(T (^(x',t')^( x ^))W = iS F (x';x\m) + iS l3 F '^(x';x\m) 



where we have used the only non-zero bilinear thermal averages 



(2.13) 



(2.14) 



(2.15) 



In principle we do not have to restrict ourselves to thermal distributions as given by 
Eq.( |2.8| ). In fact, we can allow for any such one-particle distribution function fp~ (a;) and 
the definition Eq.flgTg). 



3 External uniform and static magnetic field 



For the convenience of the reader, we summarize some of the relevant expressions in 
the case of a constant magnetic field B parallel to the z— direction in the gauge = 
(0, 0, — Bx, 0). Using k as a collective index for (n, k y , k z ), where n = 0, 1, 2, ... ; k y , k z are 
continuous, and the 7 matrices in the chiral representation, we can write the solutions in 
the form 

■,(±)/ 



where 



/W , .x 1 exp[±(-iE K t+ik y y+ik z z)] {±) 

/ (E K + k z )I n , ky (x) x 



VE K + k z 



1 



v / E K + k z 



-iy/2eBnI n _ 1 . ky (x) 
-ml n . ky (x) 


f \ 

-ml n - 1;ky {x) 
-iy/2eBnI n . ky (x) 
y (E K +k z )I n ^ 1;ky (x) J 



(3-1) 



(3.2) 



(3.3) 
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^JE K -k z 



1 



In these expressions the energy E K is given by 



( -ml n .- ky {x) \ 


{-E K + k z )I n ._ ky (x) 
y W2eBnI n _ 1 .^ ky (x) J 

( tV2eBnI n .^ ky (x) \ 
-E K +k z )I n -i.-k y {x 


V -ml n -i.- k y(x) J 



E K = Jm 2 + k 2 + 2eBn 



and the I n; k v (x) functions are explicitly written 



I K kA x ) = ( — J ex P 



— -eB \ x——^- 
2 \ eB 



V2e~B ( 



x—- 



V eB 



Here i7 n is the Hermite polynomial given by Rodrigues' formula as 



H n {x) = (-l) n e 



i„2 d r ' 



dx T 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



and we define I-i- kv (x) = 0. The functions I n -k y {x) are normalized as 

dxl n;ky (x)l m]ky (x) = 5 n)m (3.9) 
if n, m > 0, so that it is easily shown that the collection of all \l/'s form a complete 



orthonormal set. The vacuum part of the propagator Eq.(2J:) is then given by (see e.g. 
Ref.H) 



00 f duo dk dk 

S F (x';x\m) ab = — tnl\?. Z ex P[- iuJ (t' - 1) + ik y (y' - y) + ik z (z' - z)} 



n=0 

X 



(2vr) ; 



uj^k 2 -w?-2eBn + W ' *»' **' 



(3.10) 



The matrix S^n; cj, fc^, entering above is given by 
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S(n;u, k y , k z ) = 



( ml n>n 

m/ n _i jn _i 

-(uj-k z )I n . n i\p2eBnl n ^ n -x 

-W2eBnI n -i tn -(w+k z ) I n -i in -i 



-(u+k z )I njr 
iy/2eBnI n -i t . 

ml n ,n 




-iV2eBnI ntn _i 
— (u — kz)I n -i, n -i 


mln-l,n-l 



where we have defined 

Similarly we find the thermal part of the fermion propagator 

S^{x']x\m) ab = / ^o^3 * exp[-iv(t? - 1) + ik y (y' - y) + ik z (z' 
n=Q J 

x2wi8((jj 2 — kl—m 2 — 2eBn)fF(u!)Sab(ri; u, k y , k z ) , 
where /p(o;) is the thermal distribution 

f F (u)^9(uj)f+(uj)+e(-oj)fp(-uj) . 
By making use of the completeness relation 

oo 

In;k y (x)I n;ky (x') =8(X-X') , 

n=0 



one can show [17] that the propagators Eq. ( |3.10 ) and Eq.( 3.13 ) reduce to the 
propagators in the limit when the magnetic field B tends to zero. 



(3.11) 
(3.12) 

-*)] 

(3.13) 

(3.14) 

(3.15) 
free-field 



4 Propagators in thermo-field dynamics 

The thermal propagators in Eqs.( |3.10| ) and ( |3.13| ) cannot be used for a perturbative ex- 
pansion in a naive way. The reason is that the 5-functions can occur on several internal 
legs with coinciding arguments and that such expressions are not well-defined. It is known 
that such problems can be avoided be means of a correctly derived real-time finite temper- 
ature formalism where one must invoke a doubling of the degrees of freedom. There are 
several formalisms for doing that and we shall use thermo field dynamics (TFD) since it is 
very easy in the operator formalism . In TFD the propagator is obtained as the expec- 
tation value of the time-ordered product in the thermal vacuum \Op) which is annihilated 
by the thermal operators (b\ K (j3), d\ K {(3)) and their tilde partners (&a k (/3), d\ K (f3)). 
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The TFD propagator can be given by a simple expression for independent harmonic 
oscillators. We have solved the Dirac equation exactly in the external field, but in the 
free propagator the interaction between the particles is neglected. Each mode is therefore 
still an independent harmonic oscillator, but with a different frequency labeled by the 
quantum numbers (n, k y , k z ) and corresponding to a definite Landau level. Thus, in the 
derivation of the propagator we can copy the usual procedure for free particles. 

The Bogoliubov transformation between the zero temperature and thermal operators 
is given by 




and 



dx, \ ( cos$[-J -sin^ \ / d x M 

idl / l sin^ cos*£ ) \ ^xM 



(4.1) 



(4.2) 



The number expectation values in Eq. fl2.15|) imply that the coefficients in the Bogoliubov 
matrices must satisfy 

smH^J = fP(E K ). (4.3) 



We use the convention that b and b anti-commute. The Bogoliubov matrices in Eqs. (|4. 1|) 
and ( [4.2|) , and the factors of i, are carefully explained in Ref . [plU . The definition of the 



fermion propagator varies slightly in the literature and we shall for definiteness follow 
Ref. pi | . Other propagators correspond to other definitions of the thermal doublet and 



they may be computed in a similar way. Since we are not doing any higher loop cal- 
culations these conventions do not matter. We compute the TFD propagator matrix 
as 

iS T F FD (x'-x\m) ab = (O p \T 



The structure of the propagator is the same as in absence of the external field except that 
we now expand in another basis corresponding to the new energy eigenvalues. We obtain 

00 r dujdk dk z 

iSp FD {x';x\m) ah = J2 — , v , 3 z exp[-iu(t' - t) + ik y (y' - y) + ik z (z' - z)] 



n=0 ' 



1 







xS ab (n;uj } k y ,k z )U F (uj)[ ]Uj.(u), (4.5) 



7 



where 



U F {u>) = | C ° S ^ M - Sin ^ } J . (4.6) 
sini?(u;) cos 



and 



sin^H = 6(u)yfti(u)-6(-u)yJfF(-u) , 

cosi%) = - f+(u) + 6(-u)yJl - fp(-u) . (4.7) 

Here, Up(cu) is the transpose of the matrix Up{u). The Sp FD (x';x\m)^ component is, 
of course, the same as the propagator in Eqs. ( |3.10[) and (|3.13|) and the other components 
are only needed in higher loop calculations. 

The derivation of the propagator in this Section can be repeated for a non-equilibrium 
distribution if only we assume certain factorization properties of the density matrix. The 
essential assumption is that there are no non-trivial multiparticle correlations so that 
everything is determined in terms of the single particle distribution. This freedom amounts 
to replacing the functions fp(E K ) with some other positive functions that describes the 
distribution. Some applications of such a formalism in absence of the external field can 
be found in Ref . |22| . 



5 The effective action 



As was shown by Schwinger a long time ago [23], an external electromagnetic field, slowly 



varying in space and time, can be treated to all orders in the external field in the weak- 
coupling limit. Here we make use of a technique similar to that of Schwinger 's in order 
to evaluate the thermodynamical partition function in a static uniform magnetic field B 
for charged fermions as well as for charged bosons. 

5.1 QED and Charged Fermions 

The generating functional of fermionic Green's functions in an external field Z[fj, 77, A^], 
formally defined by 

Z[fj, 77, Ap] = J d[ij}d[ij} exp[2 J d^xi—F^F^ + m)ij - # + $77)] , (5.1) 
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describes second-quantized electrons and positrons interacting with a classical electro- 
magnetic field expressed in terms of the vector potential A^. The expectation value of ip 
(and x/j) can formally be fixed by choosing appropriate f) (and 77), i.e. <p(x) = (ip(x)) = 
— i5/5fj(x) log Z (and similarly for ip). The equation of motion for tp(x) tells us how the 
electrons interact with the electromagnetic field which includes effects due to all virtual 
e + e~-pairs. 

The fermionic Gaussian functional integral in Eq.( |5.1|) can formally be performed with 
the result that 



Det Jp — m)] exp 



J d A x (--F^F^ + J d 4 yr](x)S F (x;y\m) V (y) 



(5.2) 



where Sf{x; y\m) is the external field vacuum propagator as given by Eq. (|3.10| ). It follows 
that (f(x) satisfies the Dirac equation in the external field, i.e. [iip — m)ip(x) = 0. 

The functional determinant Det(i jp — m) gives rise to a contribution to the effective 
Lagrangian density £ e //. Using a complete orthogonal basis to rewrite log Det as Tr log, 
the effective action can thus be written 



eff 



d A xC 



eff 



d x 



— i Tr log jp — m)] 



We now write the effective Lagrangian density as 



eff 



C +C 



1 ) 



(5.3) 



(5.4) 



where the tree level part in the case of a pure magnetic field is 

1 



(5.5) 



and £i corresponds to the functional determinant. Differentiating Eq.( f5.3D with respect 
to the fermion mass we now find the one-loop correction according to 



dm 



i trSV(^; x\m) , 



(5.6) 



where the trace now only is over spinor indices. After a straightforward calculation of 
the trace using Eq. (|3.10|) , we obtain in terms of renormalized quantities the well-known 



result p3f that 



00 ds 



]ir 2 Jo s 3 



esB coth(es-B) 



1 - -(esB) 2 



exp(— m s) 



(5.7) 
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We have here performed the standard renormalizations leaving eB invariant, i.e. 

e 2 — e 2 (l + Ce 2 ) , (5.8) 



where the divergent constant C is given by 

C =^r~2 f°-exp(-m 2 S ) . (5.9) 

We shall now find the corresponding correction S^fi = J d 4 xC^ , to the effective 
action S e ff at finite chemical potential and temperature such that 

C e ff = C +£i+£ p e f f . (5.10) 

We notice that the correction £e//> due t° the presence of thermal fermions, can be written 
in the form 

= iTrSp^ix; x\m) . (5-H) 

dm 

By performing the trace operation in Eq. fl5.1l|) , using the thermal propagator Eq. Q3.13|) , 
we obtain 

AeB 00 2 



n=0 A=l 

oo 



X 



poo 

/ dkk 2 5(u 2 - k 2 -2eB(n + A - 1) - m 2 ) , (5.12) 
Jo 



where we have integrated by parts with respect to k. We, therefore, see that jO e y J\ is 
directly related to the partition function Z(B,T,fi) of the relativistic fermion gas in the 
presence of an external magnetic field B in a sufficiently large quantization volume V, as 



given in for example Ref . p4j , according to 



r p,n logZ(73,r, (j) 

Leff ~ pv 

eB ~ 2 A 00 „ k 2 



EE/ 

n x 1 J — OO -C/X n 



(2tt) 2 ^oZ^J-oo E x , n 
X (l+exp[/5(E A , n -/i)] + l + exp[/?(£ Ain + //)]) ' (5 ' 13) 



where 



E Xn = Jk 2 + 2eB(n + X - 1) + m? . (5.14) 
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Separating the field independent part we write 



= + , (5.15) 



J eff 

where 

3/2 



= / ^(„ 2 - m 2 )/ F (a7) -m 2 ) 1 . (5.16) 



37T 2 J-oo 

We, therefore, conclude that the field independent thermal correction to the Lagrangian 
density can be identified as 

r p,v _ logZ(T,;x) _ F(T,fi) 

where Z(T, /a) is the partition function, and F(T, /i) the free energy, for an ideal e + e~-gas 
with particle energy E = \/h 2 + m 2 , i.e. 

log Z(T,li) f d 3 k 



I wr* ( log[1 + e ~" {E ~" )] + log[1 + e ~" {E+ " )] ) ' (5 - 18) 



V J (2tt) ; 

consistent with the general identification above. Using the identity 

exp(— Id) r°° dt / 1. 9 l.\ .„ 

1 ! - ' : exp \ --{xH+ -) , (5.19) 



F| Jo y/2M V 2 V * 

the following representation of C^, valid for < m, can be derived in a straightforward 
manner 

,p,n 1 [°°ds ( (3 2 l 2 2 \ coshQg l/i) 



^r = 4^E(-l) m y ^exp^-^--m 2 S j^^[e5 S coth(e5 S )-l] . 

(5.20) 

In the case \i = 0, Eq. (|5.20| ) agrees with the result obtained in Refs.0, However, it is 



not always obvious, when written in this form, to see how to extract the physical content, 
and particularly not obvious how to generalize L e f^ to > m, since then it appears to 
be divergent. In particular we notice that the high T behaviour given in Ref.]l2| is not 
correct. As explained in Appendix A it is, however, possible to show that Eq. ( |5.20| ) is 
equal to Eq. fl5.2ip given below, which is valid for all T and \i. 

In order to calculate the thermal part C^ff of the effective action in a more useful form, 
we have to be careful with the convergence and the analytical structure. Some details of 
the calculation are given in Appendix A. We get C p e f f = + £?'", where the ideal 
gas contribution in absence of the external field B, is given in Eq. ( |5 . 1 7|) , and 
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J 

du!0(uj 2 — m 2 )/i?(c<j) 

o 

duj6{uj 2 — m 2 )fp(uj) 



4vr 5/2 



ds 



<s(ui 2 —m?) 



S 



5/2 



[seB coth(se5) - 1] 



1 -m ir \j 



7r 7rn 
V4 ~~ IB 



{uo 2 - m 2 )) 



(5.21) 



The term with the sum over n, Ci'g SC , was neglected in Ref.fT4}| and we show in Section |] 
that it is essential to keep this term in order to get the correct physical result. We may 
also use the generalized (^-function to rewrite C^ sc in a different form, sometimes more 
suited for numerical calculations 



osc 



s3/2 



V2. 



7T 



(--, mod 


u 2 — m 2 








) 




2eB 





where modL4] is a shorthand notation for A modulo 1, i.e. 

mod [A] = A - hit [A] . 

An alternative way to write Eq.( |5.22| ) is 

ds 



(5.22) 



(5.23) 



osc 



oo r \ 



E 

n=0 



m 2 + 2eB(n + s) 



x [If ( V m2 + 2eB ( n + *)) + fp ( V m2 + 2e£? ( n + S ))J C(~2> s ) > ( 5 - 24 ) 

where the various Landau- level contributions are made explicit. In addition to the 
free energy has a contribution from the thermal photons, i.e. 

<"■> 

which is background field independent since there is no self-interaction among abelian 
gauge fields. 



5.2 QED and Charged Scalars 

The formalism used so far applies also to scalar QED. We give some of the corresponding 
results here for completeness. Equation ( |5.6| ) becomes in this case 



dm 2 



-iGp(x; x\m ) 



and the thermal propagator is 

(G F (x;x\m 2 ))/3 tfM ■■ 



~ f dudk y dk z , 

2^J (9 ^ 3 ^n;k y {X) 



X 



(2tt) 3 
1 



to 2 - El + ie 



-2ui5{uo 2 -E 2 n )f B {uj) 



(5.26) 



(5.27) 
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where 



and 



El = k 2 z + (2n + l)eB + m 2 , 
9(u) 9{-uj) 



+ 



It is rather straightforward to obtain the correction 



l 



167T 2 JO S 



ds 



eBs 



— exp{—m s) 



sinh(e5s) 



1 + 



(eBs) 
6 



(5.28) 
(5.29) 

(5.30) 



to the effective action in the vacuum sector. At finite chemical potential and temperature 
we similarly find the following contribution to the effective action 



^eff 



1 



6tt 2 



du9{u 2 - m 2 - eB)f B {u)(u 2 - m 2 f' 2 



Jdu9(u 2 -m 2 -eB)f B (u) JL | 



ds 



n 5/2 J s 5/2 

x 'eB\ 3/2 



-s(u} 2 — m 2 ) 



eBs 



sinh(ei?s) 



- / du9(u 2 -m 2 - eB)f B {u) JL f (^) sin ^ - ^(u 2 - m 2 - eB) 

(5.31) 

The zero temperature part £,\ was derived in Ref.j23|]. Physically this effective action is 
quite different from the fermionic one. We shall not pursue this investigation here but 
only make a few remarks. Since for charged bosons there is no sharp Fermi surface, there 
are no de Haas - van Alphen oscillations either. Furthermore, even the energy of the 
lowest Landau level [n = 0) depends on B, so that, for example, in the case of a vanishing 
chemical potential, the number density is Boltzmann suppressed for large fields. 



6 The Physical Content of C e ff 



There are several dimensionful parameters related to C e ff, i>e. T, /i, m, and B, that can 
be large or small compared to each other. We shall discuss some of these limits which 
we think are particularly interesting. A central feature of a fermion gas is whether it is 
degenerate or not, i.e. whether or not the Fermi surface is sharp on the scale of the Fermi 
energy. With an external magnetic field it is also important to compare the smoothness 
of the Fermi surface with the spacing of the Landau levels. A criteria for the de Haas 
- van Alphen effect is that the distance between the Landau levels close to the Fermi 
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surface is considerably larger than the diffuseness or fluctuations in the Fermi surface due 
to finite temperature, electron - electron interactions, impurities etc. This can sometimes 
be achieved even at high T by having large p and B. 

The effective Lagrangian is here given as a function of the chemical potential p. In 
many situations it is more natural to consider the expectation value of charge density Q/V 
as given, where Q is the total conserved charge. It is calculated from Q/V = —ep(p), 
where 

which in the case of vanishing magnetic field and temperature reduces to 

and p has the same sign as p. For large B field this relation gets substantial correction, 
see e.g. Section |6"72| . We notice that p is equal to the difference between the electron 
and positron number densities, that may be useful on comparison with condensed matter 
physics calculations. In other situations one may consider adiabatic changes of B, and 
then keep the entropy fixed, or the pressure. All these different cases are described by 
suitable Legendre transformations of the thermodynamical potential F. 



6.1 The de Haas - van Alphen effect 

At low temperature one may attempt an expansion in T using Sommerfeld's method 



We assume that p > m since for \p\ < m the thermal contribution is exponentially 



suppressed. The Sommerfeld expansion for a function H(u) is 

d h II 

n ! 

where 



POO POO 00 

/ du f£(uj)H(Lu) = /+(m) / du H(lu) + £ T n a n 



duo 



n-l 



■i 



(6.3) 




(6.4) 

but the odd powers of T are exponentially suppressed. This formula can be applied to Cq'^ 
and Ci'r e g, but in Ci'£ sc performing the derivative inside the summation sign is not allowed 
since the sum is not uniformly convergent, and when acting on the form containing the 
^-function there will obviously be divergences at discrete points. This indicates that an 
expansion in mT/eB is not possible. Anyway, the T = part of Cx'g SC can be calculated, 
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and if we in particular assume {T = 0, eB <C /i 2 — m 2 <C m 2 } we get 



(e£) 2 V/U 2 - m 2 (e5) 5 / 2 ~ 1 



12vr 2 



m 



/ 7r u 2 — m 2 ~ 

cos mr 

U e£ 



1 

71 



(6.5) 



This is a non-relativistic limit (in the sense that the kinetic energy is much smaller than 
m) with a degenerate Fermi sea and a weak external field. 
The vacuum correction is in this limit given by 



so that the finite density correction 



{eB) 2 feB\ 



3607T 2 \m 2 ) 



therefore dominates over Ci when 



{eB) 2 /37r 2 p 



1/3 



12tt 2 



77T 



m- 



1/3 



> 



eB\ 



30(3tt 2 ) 1 /3 \ m 2j 



or equivalently, in terms of the chemical potential 



v> 2 - m 2 ^ 1 feB^ 



m 



30 \m 2 / 



(6.6) 



(6.7) 



(6.S 



(6.9) 



This is always satisfied in the limit {eB <C /i 2 — m 2 <C m 2 }. 

Even though the S 2 dominates over B 5 ^ 2 for small B, the magnetization of the heat 
bath@ gets a larger contribution from Cio SC , 



M = M rea + M osc = — 

reg osc VdB 



1 OF dClf f 



dB 



(6.10) 



where to the lowest order in the magnetic field 



ey eB{ji 



2_ m 2\ oo 



47r 3 m 



E 



/ e \ 
V2m 



n 3/2 
2 ™2 



7T 



n=l 

2 (^ 2 -m 2 )C 



/ 7r /i 2 — m 2 ' 

sin n7r 

U eB 



mod 



/i 2 — m 2 
2eB 



and 



(^) " ^ • 



re9 \2mJ 3 



(6.11) 



(6.12) 



4 The vacuum contribution to the magnetization is not included in Eq.( 6.1(j| ) since it is very small for 
small B. 
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The ^-function has its maximal modulus at C( — f j l) = C( — §> o) ~ —0.208, which implies 
that the peak magnetization from the oscillating term is larger than that from the regular 
term for {eB < 0.78(/z 2 — m 2 )}, i.e. when the approximations used here are valid. Defining 
the magnetic susceptibility as the response in the magnetization due to a magnetic field, 
i.e. M = xB, as in Ref . ||26|| , we get exact agreement with this reference, but not with 
which have an extra factor (— l) n in the sum over n, that we find only should 



Refs.|27|, 

be present in the case of spinless bosons. In Section |T2] we give an argument why our 
result has to be correct. 

The oscillatory behaviour as a function of B is well-known as the de Haas - van 
Alphen effect. The frequency of this periodic function agrees with the one derived by 
Onsager |29 ]. Equation (|5.21| ) describes the full relativistic generalization of this effect, 
and in Section [7| we consider some astrophysical applications where the non-relativistic 
approximation is not valid. The distance between the magnetic field of two adjacent 
minima of the magnetization is determined by 



1 



eB, 



2tt 
~A 



eB i+ i 

where A is the area of an extremal cross section of the Fermi sea. 

Sometimes (e.g. in Ref. ||25||) the magnetic susceptibility is defined by 

dM 

X 



(6.13) 



dB 



5.14) 



but again we find that the sum over n does not converge, and that the form containing 
the ^-function contains divergences at discrete values of B, and is poorly illuminating. 



6.2 Strong B-field 



In the limit of strong field, {eB 3> T 2 ,m 2 , |/i 2 — m 2 |}, we can see from Eq. fl5.13|) that 



only the lowest Landau level contribute and L e ft goes like a linear function of eB. We 



shall now reproduce this result from Eq. (|5.21|) and it turns out to be rather non-trivial. 
The leading B dependence in the first term in Eq. (|5.21|) is obtained by scaling out eB 
and taking eB — > oo in the remainder. The total contribution is, apart from the thermal 
integration (see Appendix B) 



3/2 



(gg) 

47r 5/2 



X 



dx 

5/2 



' x coth x 




(6.15) 
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but this is actually identically zero. The next subleading term can be shown to be 

a pR roc 

C^ = — du9(u 2 - m 2 )f F (Mu 2 - m 2 , (6.16) 
2n 2 J-oo 



which is exactly the leading term from Eq. (|5.13 ). This calculation shows that the oscil- 



latory term in Eq. ([5.21 ) is absolutely necessary to cancel the B 3 ! 2 term and to give the 



correct linear term. Also, notice that the expression presented here for this term has to 
be correct, without the extra factor (—1)™ of Refs. f28|, |27|j , for the B 3 ! 2 terms to cancel. 

In this limit of strong magnetic field the thermal and density corrections given above 
are small compared to 

The vacuum polarization effects are dominating here, which comes quite naturally, since 
the magnetization from real thermal particles becomes saturated when all spins are 
aligned, whereas the magnetization from vacuum polarization increases like B log B. This 
has not always been recognized in the literature |24| . 



Another issue when the B field is strong compared to fi 2 — m 2 is that the relation 
between p and \i is changed. In fact, we have from Eq. (|6.16|) at T = 



eB 

P(») « TTTvV 2 ~ m2 ■ (6-18) 



The linear dependence on the Fermi momentum kp = V '/z 2 — m 2 can be understood 
from the fact that only the lowest Landau level is filled and therefore the phase space is 
essentially one-dimensional. 

6.3 Weak B-field 



In Section |6.1| we had an expression for C^'Ji in a weak (-C fi 2 — m 2 ) field but T 2 still 



smaller than eB. An expansion for B smaller than all other scales would be desirable 
but there are some subtleties involved in such an expansion. The vacuum part can be 
expanded in a naive way and we get 

A = -— £ (_) (-1)*C(2^ + 2)r(2*) . (6.19) 

k=l 

This series is not convergent but Borel summable for small eB/m 2 so we expect the first 
few terms to be a good approximation for weak fields. Expanding the integrand of C^ eg 
(see Eq. (|5.2lD ) in powers of B leads to the same problem after the s- integration. Moreover, 
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the ^-integration becomes infra-red divergent, for higher order terms. We cannot even 
expand the integrand of Cio SC in powers of B, but after repeated partial integrations with 
respect to uj we obtain 

___4 oo / „ R \ 5/2+fc 

= ISiSS C(V2 + *)(-l)^ 



x (™ , ^V^(^«) + /'») 



) 



(6.20) 



where [fc/2] is the integral part of k/2. When > m the factor with derivatives of fp(oj) 
at uj = m contains powers of m/T. These factors, combined with the B/m 2 factors, 
show that we must have {B <C m 2 ,T 2 } in order for the expansion to be valid. For 
< m these terms are exponentially suppressed at small T. We thus see that there 
is an intricate interplay between B and T in such a way that when {eB < T 2 } £f'£ c 
is smaller than Cx^ eg , as well as their derivatives. However, when {T 2 <C eB}, even 
though {e£> <C /i 2 — m 2 ,m 2 }, the 5 derivatives of £x'osc are large and show a periodic 
behaviour as shown in Section |6.1| Also the expansion from Eq. (|5.20| ) is only asymptotic. 
In view of the observations above, especially the half-integer powers of B in Eq. (|6.20 ), 



it seems unlikely that the same result can be obtained in ordinary pertubation theory 
using diagrammatic techniques. The vanishing radius of convergence for the expansion 



of Ci, and the same for Ci^ eg , also including the infra-red divergences, arise due to the 
fact that we get substantial contributions to the parameter integrals when we are outside 
the radius of convergence for the series expansion of the coth(e-Bs), i.e. for large s, for 
high order terms. We will investigate this, the non-analyticity in B, and the connection 
to ordinary perturbation theory more carefully in a future project. Some weak-field 
results can nevertheless be obtained and, for instance, the magnetic susceptibility can be 
computed in the limit {B — > 0, T <C fi 2 — m 2 }. It gets contribution only from 



regi 



b->o dB z 6n 2 \ m m J 

If we further assume that { /i 2 — m 2 m 2 } and write it in terms of the Bohr magneton 
fiB = e/2m and the density of states at the Fermi surface g(fi) = my/ Ji? — m 2 /?: 2 , we find 

2 

X = g^lfl 1 ^) = XPauli + XLandau ■ (6.22) 



It coincides with the well-known result [25[ where XPauii is the Pauli paramagnetic spin 
contribution and XLandau = —\xPauii is the Landau diamagnetic orbital contribution. 



18 



Notice that in this weak field limit the thermal corrections dominate, i.e. 

*P,n _ ( eB ) 2 VP ~ m ' 2 



12n 2 



m 



> d 



(6.23) 



where L\ is given by Eq. 

6.4 High temperatures 

At high temperatures one may find an analytical approximation in the limit {T 2 ^> m 2 ^> 
eB, n = 0}, where we have that 



LI 



(eBf 
24tt 2 



log 



T 2 



(6.24) 



and we do not agree with the high temperature and weak field limit in Ref.||12||. (We 
notice the similarity between Eq.( |6.24j ) and Ci for eB ^> m 2 in Eq.( p.l7|) .) The thermal 
contribution thus dominates over C\ as given by Eq.(|6.6|) when 



T 



m 



~> exp 









\m z J 



(6.25) 



i.e. when the approximations used here are valid. 



7 Some Astrophysical Applications 

As mentioned in the Introduction, strong magnetic fields at finite temperature and den- 
sity are situations that are frequently encountered in astrophysical contexts. We have 
investigated the possibility of some interesting behaviour mainly for white dwarfs, neu- 
tron stars and supernovae since they present the most extreme conditions while still being 
directly observable, in contrast to e.g. cosmic strings, the existence of which has yet to be 
confirmed. We can use the effective action in two ways. Either we consider the response of 
the system to a given external magnetic field H, or we study the properties of an isolated 
system with only the induced magnetic field. In the first case the free energy is given by 

F = -A(B)-/2#(B), (T.l) 
where B is determined by the mean field equation 

B = H + M(B) = H + ^ + ^ L . (7.2) 
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The magnetization M(B) is thus calculated in the presence of the microscopic magnetic 
field B. Note that we include both the contribution from real electrons in the heat 
bath and virtual electrons from vacuum polarization. If we consider the dynamics of the 
system without any external field we should add £ to the effective action and determine 
stationary values of the field by 

dC eff _ 3d 9Clf f _ 

-dB--- B + ^B + ^B--°> (7 ' 3) 
which, of course, is the same as putting H = in Eq. (|7.2j ). As discussed in Section ^]2] 



the vacuum contribution is dominant for large fields. Using the result from [24] we see 
that for T = m the thermal contribution saturates at about eB = 10 m 2 . At that value 
of the magnetic field, the vacuum contribution is about twice as large as the thermal and 
cannot be ignored. 



It would be most interesting if we could find astrophysical objects showing the de Haas 
- van Alphen oscillations. The magnitude of the oscillations might then be large enough 



to effectively trap the magnetic field in a local minimum satisfying Eq.(7.3). A candidate 
for such a system is a neutron star with a strong B field and a degenerate electron gas. In 
order to get de Haas - van Alphen oscillations as a function of B the spacing of Landau 
levels near the Fermi surface need to be larger than the spreading of the Fermi surface 
due to finite temperature. If the n-th Landau level is at the Fermi surface, E n = /x, then 
we require E n+ i — E n > T. For /i 2 ^> eB, which is the case for neutron stars, we get the 
condition 

eB>fxT . (7.4) 

According to Appendix C we can even get a more stringent condition in the case of large 
chemical potential 

eB > 2<n 3 f iT . (7.5) 
As a comparision, we find in the non-relativistic case , instead of Eq.( [T.4| ) that 

eB>mT, {eB , fi 2 - m 2 < m 2 } . (7.6) 

In order to see any oscillations the field must not be so high that all fermions are in the 
lowest Landau level, i.e. integer n above must be greater than unity, that gives 

(/2 2 - m 2 )/2 > eB . (7.7) 
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White Dwarf 


Neutron Star 


Supernova 




1.02 [g 


6 • 10 2 § 


6- 10 2 § 


T/m 


2- 10~ 3 [30| 


1 1 


1 ■ 10 2 @ 


eB/m 2 


2 • 10~ 6 |, | 


2 • lO" 1 || 


2 


(/i 2 - m 2 )/(2eB) 


1 • 10 4 


2 ■ 10 6 


2 ■ 10 5 


eB/iiiT) 


1 • 10~ 3 


3- 10~ 4 


3 • 10~ 5 



Table 1: Typical values of eB,T and \x for some astrophysical objects, and an indication of 
the possibility for oscillations in the magnetization. The references are given in brackets. 



Co K) 


& (m 4 ) 


C%* (m 4 ) 


C-l'reg (^) 


C%jL K) 


-2 • 10~ 2 


2 ■ 10~ 6 


1 • 10 9 


6 ■ lO" 2 


1 • 10^ 3 



Table 2: The different parts of the effective Lagrangian for a typical neutron star, in natural 
units. 

Approximate values for eB, T and \x for what we find the most interesting astrophysical 
objects in this context, a supernova; a neutron star; and a white dwarf, are given in 
Table 1. According to above, the number in the last two rows of this table should be 
greater than unity for de Haas - van Alphen oscillations to appear. That is not the case 
in either of the situations. For a neutron star we have numerically computed the different 
parts of the effective Lagrangian, and the corresponding magnetization. The results are 
given in Table 2 and Table 3, respectively. The effective Lagrangian is totally dominated 
by the thermal contribution in absence of a magnetic field , Cq , due to the extreme 
chemical potential. We would like to stress that there are no oscillations in the so called 
oscillating part of the magnetization, M osc , in this region of parameters. Obviously we do 
not expect to see any de Haas - van Alphen oscillations unless the neutron star is very 
cold (T = 0(1) eV), or if the electron density is very low in some region, for example 
close to the surface, where the field still is strong. 



Mi (em 2 ) 


M reg (em 2 ) 


M osc (em 2 ) 


2 • 10~ 5 


4- 10~ 2 


1 ■ lO" 2 



Table 3: The different parts of the magnetization for a typical neutron star, in natural units. 
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Figure 1: The vacuum and thermal contribution to the magnetization showing de Haas - van 
Alphen oscillations as the temperature is lowered. The chemical potential is [i = Am. 



In order to investigate the behaviour of a relativistic gas of fermions showing de Haas 
- van Alphen oscillations, we have numerically calculated the effective action, and the 
magnetization for {fi/m = 4; T/m = 0.01, 0.1, 1.0}. The latter is shown in Fig. 1. 
We see clearly how the oscillations disappear as the temperature is raised. There is 
also a last oscillation at about eB ~ 7m 2 which occurs when the second Landau level 
leaves the Fermi surface. For the values above we do not find any non-trivial solution 
to Eq. (|7.3| ) because the tree level —B dominates. It is, in fact, only for a rather limited 
range of parameters that £-\'o SC can give local maxima for the total effective action. As an 
example, let us first put T = since that only enhances the oscillations. Then we look at 
small B so that the tree part is small. There is a chance that the a/B term in Eq. (|6.11|) 



can dominate. Using /i ~ m and making the approximation Q\— ~. 

V2e 3 / 2 



M osc <0.2- 



m 



0.005VB(fi 



x 



m) 



< 0.2, we get 



(7i 



2m7r z 

For this term to dominate over \M tree \ = B we need B < 10~ 5 (u, — m) 2 which complicates 
numerical calculations. Also, since the field is small the probability of tunneling through 
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the barrier between the maxima is not very suppressed and it is probably not an efficient 
way of trapping magnetic fields. At very large values of B the vacuum part eventually 
dominates over the tree level, but this is just the Landau ghost and we cannot draw any 
conclusion about any instability. 

'eff> 



Even if there are no local minima in -B 2 /2 -d- C p e f f , there may be intervals in 



B where —C\ — u^ft is concave, i.e. where the susceptibility is positive. Domains with 
different magnetization could then be formed in presence of an external field, just like in 



some solid state materials 27 



8 The Effective QED Coupling 

The charge renormalization given by Eq. ( |5.8|) also leads to the weak coupling expansion 
of the QED /3-function, i.e. 

X±a(X)=P(a(\)) = ^a 2 (\) + 0(a 3 (\)) , (8.1) 
dX 3n 

where A is a momentum scale factor. We notice that due to the scale invariance of eB, 



8.2) 



we can also define an effective coupling constant from C e ff as |23|, |3T 

1 _ 1 dC eff 

e 2 ( eJ B,/i,T) eBd(eB) ' 
that gives for the electromagnetic fine structure constant a(eB, fi,T) = e 2 (eB, //, T)/47r 

1 _ 1 1 d(d + £f ") 
a(T,fi,B) a aB OB 1 ' ' 

in analogy with the definition of the renormalized coupling in the vacuum sector in con- 
nection with Eq. (|8.1| ). Special care has to be taken when evaluating the derivative of the 
oscillating term in Eq. (|5.21|) . In the limit when eB = 0, we obtain the effective coupling 
a(T, fi) = a(T, B — 0) given by 

1 1 2 , 9(uj 2 - m 2 ) , . , . . 

TT^ = -~T d0J / 2 a frM • 8.4 
a{± , fx) a o7r J-oo y/to 2 — m 2 

When T = 0, we therefore get an effective coupling a(/i) = a(T = 0,/i) such that 



1 1 2 >o g (M + ,/i^_ 1 , (8 . 5) 



a(n) a 3i \ m V va 2 
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In the limit fi = 0, we find the following asymptotic behaviour of the corresponding 
effective coupling a(T) = a(T, // = 0), 

dx 1 12, /T\ . . 

--—log - , (8.6) 



a(T) a 3vr Jp m _ 



\2 e x + 1 a 37r \m 



for T ^> m. It is now clear that (only) for fi 3> m and T»m the effective couplings a(/i) 
and a(T) are solutions to the renormalization group equation ( |8.1| ) when A is identified 
with \i and T respectively (see in this context e.g. Refs. p2|, p~3[| ) . We also note that 
Eq.( |6.17| ) leads to an effective coupling a(B) = a(T — 0, \i — 0, B) with an asymptotic 
behaviour 



1 12, VeB\ . „, 

bg p_ , (8.7) 

that also satisfies the renormalization group equation Eq. (|8.1|) . The effective coupling 
defined in Eq.( ^.3|) can also be extracted from the residue of the thermal Debye-screened 
photon propagator (see Ref . p2| ) . 

The effective couplings as given in Eqs.(|3.5|), ( p.6|) and ( |8.7|) can be interpreted as 
follows. If we use the lowest order /3-function in Eq.( |8.1| ), then the scale dependent 
coupling is given by 

1 1 'V^ • (8-8) 



Then we can write 



a(\) a 3n \m 



1 1 2 fx\ 

-7^ log T , (8.9) 



a(x) a(X) 3tt \A, 

where x = /i, T or \feB. If A is identified with any of these scales, we can in each such 
case write 

when x ^> (m and any other scale of dimension energy). 

In terms of the effective fine-structure constant, and in the case of small chemical 
potentials, so that |//| < m, we obtain 

a{eB ,T, A = — , (8.11) 

where we have defined the functions X(eB, T, fi) = Xi(eB) + X 2 (eB, T, fi), 

dx 



1 f°° dx . wf . 
Xl (eB) = — —exp-x— 
2tt Jo x eB 



1 coth(x) + 2 

sinh 2 (x) x 3 



(8.12) 
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and 



1 00 pc 



dx 



x 



X 



exp 



(3 2 l 



2/2 



4x 



m 2 x 



coth(e73x) 



sinh {eBx) eBx 



cosh(/5//i) 



The function Xi(eS) has the following expansions 



X 1 (eB) 



eB 



, Y 

457T \m 2 J 



+ 



(m 2 ) 



if e5 <C m 2 and 



Xx( eJ B) = — tog 



m 2 



3m 2 \ 



2 eJ By 



+ 



m 
ev3 



(8.13) 



(8.14) 



(8.15) 



if eB 3> m 2 . In the case of a vanishing chemical potential we can in Eq. (|8.13| ) identify a 
$4-function, given as 



d A [z,q] = l + 2£(-l) n g n2 cos(2 



(8.16) 



n=l 



and write 



1 f°° dx . 
X 2 {eB,T) = — / — exp(-x— 
2tt Jo x eh 



x 



If ei? <C m 2 , we can write 



27T JO X 

eB(3 T 



0, exp 



4.r 



eB' 
coth(x) 

x 



sinh (x) 



4 00 / 
X 2 (ei?,T) = — ^(-1) + ^o(/9m0 - {(3ml) 2 K 2 {f3ml)0 
37r 1=1 V 



(4)' 



For T>mwe can use 



£A (x/)(-l 



--log a:; x^O , 



(8.17) 



(8.18) 



(8.19) 



to find that it leads to a log(T/m) dependence with the correct prefactor in accordance 
with Eq. (|8.9| ). (The approximation of keeping only the / = 1, as in Ref.0], excludes the 
factor 1/2, and thus is not correct.) In general we have that 

X x {eB)+X 2 {eB,T) = — / -^exp(-x— ) 

2ir Jo x eB 



-0A 



0, exp 



eB(3 2 
Ax 



coth(x) 



x 



sinh 2 (x) 



(8.20) 
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9 Discussion and final remarks 



9.1 Inclusion of interparticle interactions 

In our one-loop treatment of the effective action we have not included interactions between 
electrons. The interaction energy between the particles increases with T and /x since the 
density increases, but so does kinetic energy. For a degenerate electron gas with large 
chemical potential the kinetic energy dominates over the potential energy for electrons 
close to the Fermi surface. However, not all electrons have large kinetic energy and 
corrections from interactions have to be considered for electrons with low momenta. The 
self-energy correction for fermions at high temperature and density, but zero external 
field, has been computed in e.g. Refs.p3L p3 (in Ref.p4j only massless fermions were 



considered but it gives an indication of the correction, especially in view of the result 
in Ref.||33||). There appear some completely new collective phenomena, such as hole 



excitations ||35|| , which are not taken into account in this paper. For the particle excitations 
the dispersion relation can be approximated by an ordinary massive particle provided the 
mass is replaced by an effective T and /i dependent mass |33 



ml + 4M 2 + m t 
2 



m p = , (9.1) 



where M is the thermally induced mass which in the case of QED is 



e 2 u 2 

M 2 = ; T = 0, y. ± 

e 2 T 2 



M = ; T ? 0, fi = , (9.2) 

at least if T 3> m and /j>m. The hole excitation has a more peculiar dispersion relation 
but its spectral weight is on the other hand lower. It is difficult to make any quantitative 
estimates of the importance of self-energy corrections. We do not, however, expect that 
phenomena like the de Haas - van Alphen oscillations should be altered since it depends 
on the electrons at the Fermi surface. 



9.2 Further developments 

There are some extensions of our work that may be of physical importance. First, we 
can consider the self-energy correction of an electron in presence of an external B field. 
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From that the anomalous magnetic moment can be extracted and compared with previous 
calculations for small B field, where there appears some problems of analyticity in the 
external photon momentum at finite density |]I5| . The self-energy is also important for 



the higher loop corrections of the effective action as discussed above. The QED radiative 
corrections could effect the electroweak transition rates, relevant for the Big-Bang primor- 



dial nucleosynthesis f36fl . The photon polarization tensor should also be calculated, and 
in particular its imaginary part which is related to the decay into an e + e~-pair. Also the 
three-photon vertex is interesting since it does not exist in absence of the external field. 
Such photon splitting processes have been considered earlier in vacuum |37|, |3^, |39], [|(| 



and it would be interesting to study the correction from a thermal environment. 

The physically more complicated case of a constant (or slowly varying) E field is 
equally interesting. A plasma does not stay in equilibrium since the E field gets screened 
and the physical picture is very different from the one discussed in this paper. Yet another 
generalization would be to expand a non-constant field in powers of the derivative. Such 
an expansion has been studied in Ref.[[|l] at zero temperature. 



9.3 Conclusion 

The main objective of this paper has been to establish the correct form of the one- 
loop QED effective action at finite temperature and density to all orders in a constant 
external magnetic field, and the result differs from earlier attempts. From the form of 
presented in Eq.( |5.21| ) we have checked several limits that can be understood from 
a physical point of view. A great advantage with our expression for is that the 
thermal distribution function /f(^) occurs explicitly. This means that it is easy to study 
other thermal situations by simply replacing /f(^) with some other (non-equilibrium) 
distribution (see e.g. Ref.|22|). The importance of the thermal correction depends on the 



value of B, T and fi. In some physically interesting cases they may be large compared to 
m but often of the same order of magnitude, which makes it difficult to obtain analytical 
approximations. It is, however, possible to use Eq.( p.21|) , or the expressions in Appendix 
C, for numerical calculations. 

Even though the correction to the free energy may be small compared to the value 
without the external field there are other quantities that are effected by the presence of the 
heat bath. For instance, the magnetization of a degenerate Fermi sea shows the de Haas - 
van Alphen effect. We found, however, that for a neutron star this effect does not show up 
in spite of the extreme degeneracy and magnetic field. The reason is the relativistic form 



27 



of the energy spectrum which suppresses the oscillations at a large chemical potential. 
We also briefly discussed the importance of including the vacuum contribution to the 
magnetization when the B field is comparable to m 2 /e. 

We have, furthermore, calculated an effective coupling constant defined from the 
derivative of Cftf with respect to B. It satisfies asymptotically a naive zero tempera- 
ture renormalization group equation where the renormalization scale is replaced by T, \i 
or yfeB. 
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Appendix A 



In this appendix we give some details of how to calculate the effective action in Eq.( |5.21 ) 



First we show that Eq.( 5.20| ) is equal to Eq.( |5.21 ). To do that we start with a Poisson 
resummation in I using 

oo j2 co i 

£(_!)' exp(- -) = v / 4^^exp(-a7r 2 (2/ + l) 2 ) - - , (A.l) 
1=1 4a i=o A 

and rewrite the sum over I as a contour integral by means of the formula 

E/&2n + Q) = -g-/ • (A.2) 

The integration contour C is chosen to go from oo + it to e in the upper half plane and 
back to oo — it in the lower half plane (i.e. to G {oo + it — > e — > oo — it}), without 
encircling the origin. In this way all the poles on the positive real axis are encircled. 

We would now like to deform the w-integral to the imaginary axis and the s-integral to 
the negative real axis. This is not straightforward since there are poles on the imaginary 
s-axis and the section at infinity has to bo chosen to give a vanishing contribution. It is, 
therefore, necessary to divide the integral into several pieces and to do the deformation 
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for each piece separately. Let us start with the part where u is in the upper half plane. 
Then the s-contour can be deformed to the negative imaginary axis, but to the right of 
the poles. After that the u contour is deformed to the positive imaginary axis. Finally, 
for \u\ > m we further continue the s-integral to the negative real axis and pick up the 
poles on the negative imaginary axis, while for |u;| < m we deform the s-contour back to 
the positive real axis. 

The whole procedure can be repeated for u below the real axis, reflecting all deforma- 
tions around the real axis. To get the correct convergence for the cu-contour deformation, 
the constant —1/2 in Eq. flA.l ) should be associated with the u in the lower half plane. 
After summing the pieces there is only a contribution from \u\ > m, as expected, and it 
consists of an s-integral and a sum over the residues of the poles. 

In the deformations above we have been careful with the convergence for large |s| and 
|c<j|, but we have said nothing about the possible singularity at s = 0. One way of dealing 
with that is to multiply the expression with s v and perform the integration for such a v 
that there is no divergence at s = 0, and to do the analytic continuation at the end. 

Equation ( J5.21 ) can also be obtained from the thermal propagator in Eq. ( |3.13[ ) by 
representing the 5-function as 

2mS(x) = ilm = ilm / ds e~ l{ - x - le) . (A.3) 

x — ie Jo 

Then the k y and k z integrations can be carried out (using Eq. (|3.9|) as well). The summa- 
tion over n is just a geometric series but it is not absolutely convergent so we sum only 
to a finite A" and take the limit A" — > oo at the end. This gives 

N->oc 7T 6 ' Z J_oo 27T JO S 1/Z ~ u 



X 



-oo 

A2sB nA2NsB 



1 + e l2sB 2e l 



(A.4) 



J _ e i2sB ^ _ e i2sB 

where we also have introduced the dimensional regularization v in 4 — 2v dimensions, 
and we are to analytically continue to v = in the end. Keeping v large enough that 
the integral is absolutely convergent, the expression above can easily be integrated with 
respect to m to yield C? e ft- To be more precise, there is an integration constant from the 
lower limit in 

% I dm'Tr Sp^(x; x\m) = C^fAra) — £fj-j(m ) . (A.5) 

J mo 

We expect that in the limit m — ► oo the thermal part of the effective action is zero since 
an infinitely massive particle has zero Boltzmann weight. Therefore we let m — > oo and 
thereby put the integration constant C^(m Q ) to zero. 
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The poles in the last factor in Eq.( |A.4j ) cancel for finite N, and we cannot let N — ► oo 
in a naive way before deforming the s integration contour to the imaginary axis. The two 
terms have to be treated separately so we must choose an integration contour for s slightly 
above or below the real axis. Since, according to the discussion above, C^(m — > oo) = 
we see that the the original contour must be chosen slightly above the real axis. Depending 
on the sign of uj 2 — m 2 (or uj 2 — m 2 — 2eB(N — 1) in the the second term) we deform the 
s-contour to either the positive or negative imaginary axis. In one of the cases we get a 
contribution from the poles. After deforming the contours we take the N oo limit and 
also take the limit v — > what concerns taking the imaginary part, in order to get a more 
apparent expression, but we still need to keep v > to have the integration over s finite, 
with the result 

1 r°° ds 

Tin 



^eff 



dujQ(ui 2 — m 2 )f F (uj) 
du8(uj 2 - m 2 )f F (uj) 



n=l 



. ,^ v e- s{u> ~ m ] seB coth(seB) 

3/2 



n J 



M111 4 



eB 



[uj 2 — m 2 ' 



. (A.6) 



Actually we must have v > 3/2, i.e. less then one dimension, but we may just consider it 
as an analytical continuation in u, in order to be able to change the order of integration. 
If we now take the limit B — > 0, we get 



1 



4vr 5/2 
1 

AtTJ 2 



(Lj9(uj 2 
dujQ[uj 2 



m 2 )f F {uj) / ds s -5/2 e -s(- 2 - m 2 ) 



m 2 )f F {uj){uj 2 - m 2 fl 2 - v Y[v - 3/2) 



(A.7) 



We may now take the limit v — > to get Eq.( |5.17] ), and after subtraction of this term we 
may also let v vanish in Eq. ( |A.6| ) and get Eq.( |5.21| ). 

Appendix B 

In the limit of very strong fields {eB 3> T 2 , m 2 , \fi 2 — m 2 |}, the first term in Eq.( |5.21| ) can 
be written as 

<V% = £^ 2 --^)^^/ o °°^(^oth S -l) . (B.l) 
Similarly we find in this limit 



C(% c = - / duj6(uj 2 - m 2 )f F (uj 

J — oo 



(eB) 



3/2 



■am 



(B.2) 
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where ( is the Riemann ^-function. It can be shown by residue calculations that 

r^ (scoths - i)= /! c(3/2) ' (R3) 

so that the 0(B 3 ^ 2 ) terms cancel in this limit. In order to extract the next term in the 
strong field expansion of we consider the expression entering in the to integral in 
Cf'^, expanded for large B, i.e. 

w { f ~ m2)l " 1)(s coths " 11 

2 JS, 1 / r 7r nn , 9 7T\] . 

-^E^Hi-^" 2 -» 2 )]— i)} ■ < B - 4 > 

If we now use the cancellations depicted above, and the fact that the sum converges 
towards an integral in the limit B — > oo, the expression above may be written as 




By performing the integrals in this expression, we find the following leading contribution 
£f " = ^- / ^(cu 2 - m 2 )/ F (o;)Va; 2 - m 2 . (B.6) 

Z7T J — oo 
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Appendix C 



In the case of large chemical potentials, e.g. in a neutron star, when (/i 2 — m 2 )/2eB 3> 1, 
the form for £i'o SC given in Eq.( p.22[ ) is difficult to handle due to the rapid oscillations in 
the ^-function. Let us instead start from Eq.( |5.21| ), and rewrite it as 



osc 



ni 



4 {eB\ z ' 2 ™ 



27r 3 \m 2 J 
where we have defined 

L = 



m < exp 



71=1 



7r irn 
v 4 eB jvn? 



dx- 



ex P 



(C.l) 



(C.2) 



1 + exp[m/3(x — /V m )] 
Since the exponential function here is oscillating rapidly and we desire a rapidly decreasing 
function instead, we close the contour with a circular section at infinity, a straight line 
from the origin to infinity with complex argument 7r/4, and the small section from the 
origin to x = 1, and use Cauchy's theorem to get 



3 27r/4 



dx 



exp 



irn 2 



eB/m 



o 1 + exp[m/3(e J?r / 4 a; — fi/m)] Jo 1 + exp[m[3(x — fi/m)] 



dx- 



exp f 



- nn 2 



eB jm? 



I jpole 



(C.3) 



The contribution from the residues at the poles is 

e5 



jpole 



T 

-2m — exp 

m 



-2tt 2 ^-— + mn±-— 



eB 



ex p 

J y=0 



-2n 2 n—2is - ivr 3 — (2i/ + l) 2 
e£ eB y J 



(C.4) 



where we have defined v max as the number of poles encircled by the contour 

v max = int ^ - ^ . (C.5) 

In the case of large chemical potential compared to the temperature and the square root 
of the magnetic field, we may assume the thermal distributions to be unity, and perform 
the integrals with the result 



>/4 L 



I eB/m 2 



n 



1 — erf 



n 



eB/m 2 



3 — i7r/4 



+ P n oles + O[e 



(C.6) 



It turns out that the phase from one minus the error function in Eq. (|C.6|) cancels the phase 



from exp —i (~ + 



eB/m? 



in Eq.( |C.l| ), when taking the imaginary part. The oscillations 
are thus only originating from the residues at the poles, that all have Ke[u} = /i, i.e. 
they are lying at the Fermi surface. Also, notice that the contribution from these poles is 



exponentially suppressed as exp 



-27T 2 4 

eB 



in agreement with the general discussion on 



de Haas - van Alphen oscillations in Section [7] 
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